arXiv:1507.07085vl [gr-qc] 25 Jul 2015 


Cosmological-Billiards Groups and self-adjoint BKL Transfer Operators 

Orchidea Maria Lecian^’0 

^ Sapienza University of Rome, Physics Department and ICRA, Piazzale Aldo Moro, 5 - 00185 Rome, Raly 

Cosmological billiards arise as a map of the solution of the Einstein equations, when the most 
general symmetry for the metric tensor is hypothesized, and points are considered as spatially de¬ 
coupled in the asymptotic limit towards the cosmological singularity, according to the BKL (Belinski 
Khalatnikov Lifshitz) paradigm. In 4 = 3 -I- 1 dimensions, two kinds of cosmological billiards are 
considered: the so-called ’big billiard’ which accounts for pure gravity, and the ’small billiard’, which 
is a symmetry-reduced version of the previous one, and is obtained when the ’symmetry walls’ are 
considered. 

The solution of Einstein held equations is this way mapped to the (discrete) Poincare map of a 
billiard ball on the sides of a triangular billiard table, in the Upper Poincare Half Plane (UPHP). 

The billiard modular group is the scheme within which the dynamics of classical chaotic systems on 
surfaces of constant negative curvature is analyzed. The periodic orbits of the two kinds of billiards 
are classihed, according to the different symmetry-quotienting mechanisms. 

The differences with the description implied by the billiard modular group are investigated and 
outlined. 

In the quantum regime, the eigenvalues (i.e. the sign that wavefunctions acquire according to quan¬ 
tum BKL maps) for periodic phenomena of the BKL maps on the Maass wavefunctions are classified. 

The complete spectrum of the semiclassical operators which act as BKL map for periodic orbits is 
obtained. 

Differently form the case of the modular group, here it is shown that the semiclassical transfer 
operator for Cosmological Billiards is not only the adjoint operator of the one acting on the Maass 
waveforms, but that the two operators are the same self-adjoint operator, thus outlining a different 
approach to the Langlands Jaquet correspondance. 

PACS numbers: 98.80.Jk Mathematical and relativistic aspects of cosmology- 05.45.-a Nonlinear dynamics 
and chaos 
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I. INTRODUCTION 

The solution to the Einstein field equations in the asymptotic limit towards the cosmological singularity corre¬ 
sponds, within the BKL (named after Belinski, Khalatnikov and Lifshitz) Ref. 1, Ref. 2, Ref.3, Ref. 4, Ref.6 
Ref. 7 paradigm, for which space gradients are neglected with respect to time derivatives, under the most general 
assumptions for the symmetries of the metric tensor, to the asymptotic limit towards the cosmological singularity of a 
Bianchi IX Universe. This way, the Einstein field equations for a Generic Universe correspond to a system of ordinary 
differential equations, whose symmetry for homogeneous universes and for inhomogeneous universe, T 2 {PGL( 2 ,C)) 
and c, respectively, allow one to define the corresponding billiard systems on the Upper Poincare Half Plane, and the 
corresponidng billiard maps for the discretized dynamics Ref.ll, Ref. 12, Ref. 13. 

The conjugacy subclasses of the Billiard Modular Group define the T2{PGL{2, G)) congruence subgroup of PGL{2, G). 
The discretization of the dynamics allows for the definition of language codes, i.e. the composition of transformations 
that describe the continuous billiard dynamics and according to the time evolution of the Einstein field equations; 
the analyses in therefore cannot be applied to Gosmological Billiards. 

The unquotiented cosmological billiard groups therefore describe the (time) evolution of the dynamics, for which 
the insertion of any composition of transformations corresponding to the Identity within the unquotiented maps 
would correspond to (with respect to the time evolution and therefore with respect to the symmetries of the solution 
of the Einstein field equations) unphysical sequences of trajectories. 

The eigenvalues for the BKL quantum operators that constitute the quantum maps and the semiclassical ones 
are defined according to the language codes of the big billiard and of the small billiard. 

Cyclic identities for perioidc orbits of the big billiard map and of the small billiard map define the parity of the (with 
respect to the corresponding WDW equation) suitable Maass wavefunctions and therefore exaclty solve the Selberg 
trace formula for cosmological billiards. 

The paper is organized as follows. 

In Section ini the Billiard Modular Group is defined. 

In Section uni the language code and the definition of periodic orbits for the Big Billiard Group are defined. 

In Section HVl the language code and the definition of periodic orbits for the Small Billiard Group are defined. 

In Section |Vl a comparison of the systems is outlined. 

In Section rvTl the quantum regime and the semiclassical transition are analyzed by the definition of the semiclassical 
Poincare surface of section for cosmological billiards and the implementation of the Selberg trace formula according 
to the sign acquired by the eigenvalues of the quantum BKL maps for periodic Cosmological Billiards orbits. 
Outlook and Perspectives arebirefly stated in Section IVlII 
Brief concluding remarks follow. 


II. THE BILLIARD MODULAR GROUP 

The Billiard Modular Group (BMG) is defined on the asymmetric domain delimited by the geodesics A, B, C, such 
that 


M : u = 0, (2.1a) 

B:u=l (2.1b) 

C : = I, (2.1c) 

for which the following transformations are defined 

Az = —z, (2.2a) 

Bz = l-z, (2.2b) 

C = i. (2.2c) 


a. The language eode for the BMG According to these transformation for the asymmetric domain (12.IL any 
matrix of the BMG can be written as one of the following 


• A, B,C, AC, BC, [BCf-, 




• I, Sr, T’^A, BT^A, BA; 
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where the matrix I is the identity matrix, the matrix T is defined as T = AB, such that Tz = z — 1, T° is its n-th 
iteration, the matrix M can be one of the following: /, T°, ,BT“, CB, B; the matrix N can be one of the following: /, 
C, A, ^T“, CA. This way, the language code for the BMG is defined. 

By a suitable transformation (conjugation), any matrix of the BMG can be written as one of the following 

• elliptic transformations A, B, C, AC, BC, {BC)"^; 

• parabolic transformations T°; 

• hyperbolic transformations CT“‘CT°^..CT“'‘. 

b. Action of the BMG on the oriented endpoints The periodic orbits of the MBG can be described by hyperbolic 
transformations by a sequence of integers ni,n 2 , ■■■, rife, after identification of circular permutations, for which reduced 
matrices are defined. 

A periodic orbit of the BMG is defined by two oriented endpoints of a geodesics, — 1 < x' < 0 and x > 1, which are 
invariant under the action of hyperbolic transformations, which define a quadratic equation with integer coefficients, 

i.e. 


CT'iiCT“=...CTi"‘x' = x', (2.3a) 

CT'iiCTn=...CTi"‘x = X, (2.3b) 

such that the continued-fraction decomposition of x reads 

X = ni H- ^- , (2.4) 

n2 + -j 


the continued-fraction decomposition of x' reads 

— x' = ni -\ - ^ -. (2-5) 

”2+ ...nfc+x' 

The code of matrices, the continued-fraction expansion of the oriented endpoints of the geodesics and the expression 
of the roots of the quadratic forms are equivalent for the definition of the conjugacy subclasses of the BMG, and can 
be used to gain information one from the others. 


A. Comparison with the SL(2,Z) group 


The group SL(2, Z) is defined on the symmetric fundamental domain delimited by the sides Ai, A 2 , A 3 , described 

by 


with the transformations 


The sides are identihed as 


Ai : M = —i, ( 2 . 6 a) 

A 2 : M = 5 , ( 2 . 6 b) 

A 3 : u^ + v^ = l, ( 2 . 6 c) 

T{z) = z + 1, (2.7a) 

S{z) = -i (2.7b) 

T : Ai —>■ A 2 , (2.8a) 

S : A 3 —y A 3 . (2.8b) 


In particular, it is straightforward verified that the transformation S in Eq. (I2.7bl) acts in Eq. (j2.8bl) by identifying 
the u > 0 part of the side A 3 in Eq. (I2.6cl) with the m < 0 part of the same side, and viceversa. 
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The asymmetric domain of the BMG (12.1|) is obtained by a suitable desymmetrizzation of the symmetric domain of 
the group SL{2, Z). 

By comparison with the transformations that define the BMG in Eq. (1^ . one learns that the action of the trans¬ 
formation C in Eq. (I2.2cl) can be interpreted as a symmetry-quotienting mechanism induced on the desymmetrized 
domain m of the BMG. This symmetry-quotienting mechanism defines the language code of the BGM, such that 
periodic orbits of the BMG are classified according to this convention. 


III. THE BIG BILLIARD GROUP 


The big-billiard group (BBG) is obtained from the big billiard table, i.e. a domain defined by the three sides a, b, 

c, 


a : u = 0 (3.1a) 

b : u = —1 (3.1b) 

c : + u + = 0, (3.1c) 

for which bounces against the billiard sides are expressed by the following transformations on the UPHP 

Az = —z, (3.2a) 

Bz = -z-2, (3.2b) 

Cz = -^. (3.2c) 


The unquotiented big-billiard map 7” consists of a suitable composition of the transformations (13.2^ : the sequence of 
this composition is obtained by the continued-fraction decomposition of the variable. 

Periodic orbits of the big billiard are a phenomenon which is more complicated than its symmetry-quotiented versions 
: there exists a particular Kasner transformation fc* (which depends on m and on the considered periodic orbit of 
T) such that 

T™(m~,m'*') = m"*") . (3.3) 

The set of six Kasner transformations is a realization of the Sy, permutation group (of order 3! = 6). In fact, this 
permutation group consists of the identity, 3 transpositions [(12), (23) and (31)], and 2 cyclic transformations [(213) 
and (321)].We recall that the order of a particular group element, such as /c*, is the smallest integer p such that 
= ky. As a transposition is of order 2, and a cyclic permutation, (123) or (321), of order 3, we see that the order 
p of /c* must be equal to p = 1, 2 or 3. Therefore, by iterating (13.311 . we get 

T^^u-,u+) = kliu-,u+) = {u-,u+), (3.4) 

and mp will be the smallest such integer. In other words, is the initial point of a periodic orbit under the 

unquotiented billiard map T, with period pm, where p = 1, 2, 3 is the order of /c*. 


A. Cyclic identities 

Given any x, y in A, B, C, periodic orbits are defined by imposing the condition (as alternatively solved in (4^ 1 

T(X,y,ni)z = z 

rii 

where the epoch map Tf^x, y; rii) for a rii-epoch BKL era of the xy type is defined as 

rij — 1 

• T(^x, y, n,) = y{xy) 2 , m odd, 

rij — 1 

• T(^x, y, n{) = y{xy) 2 x, m even, 
with rii = xnp = n. 


(3.5) 
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IV. THE SMALL-BILLIARD GROUP 

The small billiard is delimited by the sides G, B, R, defined as 


G: M = 0, 

(4.1a) 

B : u = —i, 

(4.1b) 

R: u^+v^ = l. 

(4.1c) 


The transformation that describe the bounces of the billiard ball against the sides of the small billiard table (14.11) are 
those that define the SL{2, Z) group, i.e. 


i?i(z) = -z, 
R2{z) = -Z+1, 

Rsiz) = i, 


(4.2a) 

(4.2b) 

(4.2c) 


where no identification among the sides is present, and no symmetry-quotienting mechanisms for the R side of the 
small billiard table dSI) is assumed. 

The action of the transformations (14.21) on the oriented endpoints of the oriented geodesics that define the trajectories 
of the billiard ball on the small billiard table is obtained by imposing u = 0 on the UPHP variable z = u + iv, and 
results in the small-billiard billiard map t, which acts diagonally on the reduced phase space variables u = as 


Gu = —u, 

Bu = —u — 1 
Ru = i. 


(4.3a) 

(4.3b) 

(4.3c) 


Epochs on the small billiard table are defined as any trajectory joining any two walls of the small billiard sides; eras 
for the small billiard are defined as a succession of epochs starting from the side R. 

c. The small-billiard map The CB-LKSKS map for the small billiard, tcB-LKSKS^ is defined by two different 
kinds of transformations, i.e.. 


t^’‘^z = T ^SRiT "+^z, for(u''',u ) G S'j^and(u+,u ) G 5^^, 

t^'’3’3'z = T-^SRiT-^+^R^z, for(u+,u“) G (w+,it“) G 5^3^, and(w+, w") G 5, 

which act on the subregions of the reduced phase space Sl^ and defined as 


(4.4a) 

(4.4b) 


S'! • 

^ba ■ 

U < —> Ua(M+), —$ <U < —1, 

> M.^(u+); 

C2 . 

^ba ■ 

— <I> < U~ < —1, Uaiu^) <«■*■< 


St-- 

u~ < —2, 0 < «■*■ < Ua{u^), —2 < U~ < 

— <«■*■< Ua{u'^); 

^ba ■ 

— 2 < u~ < —< M+ < ?i^(u+). 

— $ < U~ < —1, Ua{u^) < 

^ba ■ 

2 < u~ < —1, 0 < < it/ 3 (it+); 



/ + ')• 


where the functions 
/+ 


• Uol[U^ 




• up{u+) : u+ = 

• u^{u+) : u+ = 

are defined. In particular, the function corresponds to the image of the function Uaiw^) according to the 

transformation B. 
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d. The language eode for the small billiard group The following language code for the small billiard is obtained 

1. R, B, G, BR, BG, RG, RB, GR, GB; 

2. RBR, BGR, RGB; 

3. Txf- 

As one can straightforward verify, the sequence RGR is not allowed. 

The trajectories Txv" are a succession of epochs such that the first epoch starts form a R wall, and the last epoch 
ends on a i? wall; between these two epochs, bounces between the G and the B wall take place, such that n trajectories 
are present. 

As classified in [], the reduced phase space for the small billiard table is characterized curvilinear domains. In 
particular, is is possible to further divide these domain according to the preimage of the transformations (14.31) on the 
RG subdomain and on the RB with respect to the BR subdomain and to the GR one. 

This way, the trajectories Txv" are defined as a succession of transformations of the kind Txf" = RXY...X'Y', 
where the transformations X, Y, X', Y' can be B or G. More in detail, they are defined on the small billiard reduced 
phase space as 

1- T = RGB...BG, n odd, u'^ > l,u~ < —1/2, {ua{u'^),ui 3 {u'^), 

2. T jiGr, = RGB...GB, n even, > 1, < —1/2, (ua(u“'"), Ui 3 {u'^), )(«“'"); 

3. Tjjsn = RBG...GB, n odd, u+ < 0, W > -1/2, {uaiu+),u/ 3 {u+),Uf'{u+), UJf{u+)); 

4. T= RBG...BG, n even, < 0, > —1/2, (ua(M+), C/”+^)(u+); 

The functions uj/(u+), Uff{u'^) are defined as 

m ■ ■)/+ — 1 -2nu~+2u~+n^-2n+5 . 

1- “ —2 -2u-+n-l ’ 

a + i ■ - 7 /+ 1 —4n+7—2nu~+4M~+n^ . 

J ■ a 2 _2_2ii-+ra ’ 

a Tl^(n+') • 11 + — 1 2nu~-2u~+n'^-2n+b . 

Uj^ya J . a — ^ 2u-+n-l ’ 

• t/”(u+) : U+ = -4 3+2^~+2»\ 
a \ / 2 2u +n ’ 

and correspond to the preimage of the RB and RG regions of the reduced phase space according to the pertinent 
combination of transformation G and B. 

e. Action of the SBG on the oriented endpoints. Periodic orbits for the SBG are defined according by imposing 
that the oriented endpoints obey the condition 

T xY^u = u. (4-5) 

According to the classification of the sequences Txv", the following quadratic equations with integer coefficients are 
obtained for the endpoints, respectively: 

1- T RG”-, n odd, — mu + 1 = 0, with to = ^ — 1 ; 

2. T RG”, n even, — mu + 1 = 0, with m = ^; 

3. T RB”, n odd, + mu +1 = 0, with to = 

4. Trb”, n even, u^ + mu — 1 = 0, with m= 

These transformations are always hyperbolic, since their discriminants reads 

1. A-rn = Vmn? + 4 for tRC” and Irb”, with n even; 

2. A-m = V m? — 4 for tRC” and Irb”, with n odd, 

as the minimum number of epochs in each T xy” era is n < 3. 

The periodic orbits defined by an even number n of epochs allow for a continued-fraction decomposition analogous 
to that obtained in the case of the Golden Ratio and of the ’silver ratios’ for the big billiard. In the case of the big 
billiard, the phase-space points that define periodic orbits identified by these ’ratios’ are placed along the function 
(«■*") of the starting box 

In the case of an odd number of epochs, this decomposition does not hold any more. Furthermore, it is not possible 
to define any transformation able to map these trajectories to this kind of decomposition. 
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/. The epoch map for the unquotiented small billiard Collecting all the ingredients together, it is possible to 
generalize the content of the sequences T xv^ and to establish a map for the unquotiented variable u = u^ relating 
each first epoch of the small-billiard eras to each last epoch of the small-billiard eras, denoted by the phase-space 
points Up = and up = u^, respectively, as 


Ml = T gy’^uf = (-1)"(m - m), 


(4.6) 


and 


Ml = T by^uf = (-1)"'(m -I- m). 


(4.7) 


with m defined in the above. 

The era transition map is obtained by composing the epoch map with the transformation S, which accounts for the 
bounce of the R, side, such that the first epoch of the successive era is defined by the phase space points of the reduced 
phase space u' = , i.e. 


u' = STxy' 


(4.8) 


V. COMPARISON WITH THE BMG 

As evident, there are several differences between the dynamics predicted by the action of the BMG, and that 
described by the SBG. 

The main difference between the BMG and the SBG is the absence of symmetry-quotienting mechanisms. In 
fact, the transformation 1^.2cl of the BMG does not coincide with the transformation ?? of the SBG. The two groups 
are obtained via two different procedures: while the BMG is obtained from a desymmetrizzation of the domain of the 
SL{2, Z) group, the BMG is due to the symmetry-quotienting of the big billiard table according to the presence of 
the symmetry walls. 

As far as the dynamical properties of the small billiard are concerned, it is important to remark that, according 
to the desymmetrized shape of its domain, no identification is possible between the sides of the small billiard table, 
such that the different periodic phenomena described by the different txv" maps cannot be identified. 

On the other hand, as far as the geometrical properties of the small billiard are concerned, as already analyzed in 
[13, the (several) geometrical transformations that would allow one to recover squared subdomains for the reduced 
phase space of the SBG do not correspond to any succession of the language code for the SBG. 

The BKL map for the small billiard consist of a different number of Weyl reflections, according to the different 
subregions of the reduced phase space where the first epoch of each small-billiard era is issued from. This property 
further explains the physical content of the Txv" maps: as a different succession of matrices is implied for every 
Txf") in particular, those with n odd will contain and odd number of (Weyl) reflections, while those with n even 
will contain an even number of (Weyl) reflections. 


A. Time evolution and cyclic identities 

The Selberg trace formula will not therefore be described according to a suitable conjugacy subclass of the modular 
group, being T2{PGL{2,C)) and PGL{2,C) larger (with respect to the modular group) groups. 

The cyclic identities of their language codes describing periodic orbits define the parity of the wavefunction solving 
the WDW equation for this implementation of the Einstein field equations, according to the number of BKL epochs 
contained in all the BKL eras which the periodic orbit consists of. 


VI. QUANTUM REGIME 

The semiclassical Poincare return map [13 , [43 , is defined as 


TEtfiq) = J^TE{q',q)'ip{q')d^q 


(6.1) 




where the integral is performed on a surface of section E, and is extended to the corresponding degrees of freedom of 
the phase space, with te a billiard map obeying the consistency equation 

0 = (1 + te)(1 - te) 

and defines the Selberg C, function according to the generalized determinants det{l — te) and det{l — te) as 

C(s) = det{l + TE)det{l — te) = 0{s)det{l — te) 

The expression of (16.11) for the case of cosmological billiards is then expressed according to the angle 9 defined in 
Figure [2] as @ , or by its expression as a function of u*, i.e. the value u = const that defines a generic Poincare 
section different from the billiard sides, and which parameterized the energy-shell reduced Liouville measure. 

The classification of all the periodic orbits according to the different maps allows one to reconstruct the complete 
spectrum of the operator te, which is the same as its self-adjoint operator Ue, which acts on the Maass waveforms. 
The density levels dE of the quantum systems and their classical periodic orbits are related at the semiclassical 


transition by the Selberg trace formula 


dE = dE + dE^^^^, 


dEosc =:J2Ar,{E)e^, 
n 


being dE and dE°‘^‘^ the average (i.e. on containing singular orbits) contribution and the one corresponidng to perioidc 
orbits, respectively, with Sn expansion of the action and An the corresponidng coefhcient within hte sum on perioidic 
orbits n. The number of epochs in each BKL eras constituting periodic orbits for cosmological billiards define the 
parity (the sign) of the eigenvalues for the quantum-mechanical description of cosmological billairds on the UPHP, 
corresponding to the BKL (towards the cosmological singularity) asymptotic limit of the WDW equation, solved by 
the (Maass) wavefunctions for cosmological billiards 


VII. OUTLOOK AND PERSPECTIVES 

g. Outlook The mathematical analysis of the discretized dynamics of Cosmological Billiards has to be considered 
as relevant for the physical characterization of the quantum regime, the semiclassical transition and the classicalized 
states of a generic Universe with respect to the present observed Universe. 

The definition of such characterization is needed for the comparison of the experimental evidence providing definition 
about the evolution of the Universe and the external (i.e. ’on the r.h.s. of the E.f.e.’s) which have to be supposed 
to have taken place as modifying the oscillatory behavior of a Generic Universe with respect to the present observed 
values of anisotropy, as well as for the anisotropy rates of the statistical distribution of matter densities as far as the 
investigation on Astrophysical scales is concerned: the introduction of isotropization mechanisms of chaotic models, 
and those of quasi-isotropization [5ll | , can allow for a comparison with observational evidence, as well as the hypothesis 
of some inflation-generating mechanism (s^ . 

The feaures of the spectrum of the energ y le vels of the wavefunction allow one to test the phenomenological effects 
obtained in Quantum-Gravitaty models [^, as far as the possible deformations of the background (geoemtrical) 
space are concerned, and allow one to compare the effectiveness of such Quantum-Gravity motivated investigations in 
modifying the chaotic properties of the billiard systems with the effects of classical scalar fields and vector fields ■ 
The introduction of isotropization mechanisms of chaotic models, and those of quasi-isotropization can allow for 
a comparison with observational evidence, as well as the hypothesis of some inflation-generating mechanism [^ . 

A complementary investigation line is constituted by and [s^. 

h. Perspectives The analysis of the solution of the iterations of Gosmological Billiard Maps which define peri¬ 
odic orbits for cosmological billiards is fundamental for the comparison of the known results for the mathematical 
properties of perioidc trajectories on domains on curved hyperbolic spaces, whose dynamics does not result as a 
symmetry-quotienting of a pre-existing dynamical system. 

The implementation of quantum statistical maps for Gosmological Billiards is based on the comparison between the 
geometrical properties of these systems, which provide the analysis with a suitable group-theoretical structure, and the 
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symmetries of the dynamics, for which the symmetries of the metric tensor define a ’smaller’ class of transformations, 
which characterizes the statistical description. This method of investigation, the Jacquet-Langlands correspondence, 
is framed in the broader Langland programmejd^ Indeed, for cosmological billiards in 4 = 3 + 1 spacetime 

dimensions, the the semiclassical Poincare return map defined by Eq. (EH) is defined by a map te for a given energy 
level corresponding to the classical configuration of energy E (at which the classical reduced phase space corresponds). 
More in detail, the operator te corresponds to any of the classical billiard map defined for cosmological billiard, i.e. 
either the big billiard unquotiented map, or the big billiard Kasner quatiented maps, such as the BKL epoch map, 
the BKL era-transition map and the CB-LKSKS map, or the small billiard unquotiented map, or the small billiard 
BKL map. 

The expression of (EH for the case of cosmological billairds is then expressed according to the angle 0 defined in Figure 
?? that defines a generic Poincare section different from the billiard sides, and that parameterizes the energy-shell 
reduced Liouville measure. 

The classification of all the periodic orbits according to the different maps allows one to reconstruct the complete 
spectrum of the operator te- 

Buondary conditions for cosmological billiards have already been thoroughly discussed in the literature. Both Neu¬ 
mann and Dirichlet boundary conditions have been proposed and motivated, according to different features that had 
to be described. 

From (??), one learns that the two different conditions, i.e. (1 -I- te) or (1 — te) correspond to Neumann boundary 
conditions and to Dirichlet boundary conditions, which correpond, on their turn, to odd wavefunctions or to even 
wavefunctions. 

It is crucial to remark that the identification of Eq. EH to a semiclassical version of the BKL map operators 
for a fixed energy shell, and the identihcation of Eq. (??) with the choice of boundary conditions is restricted to 
either the surface of section E corresponds to one side of the billiard (for the purposes of this analysis, the side b 
of the big billiard), or the surface of section E does not correpond to a side of the billiard table, but the knowledge 
of both the unqutiented dynamics and the requested maps allows one to recast the proper geodesic flow through E, 
i.e. the one corresponding to that bouncing onto a side. The classical description of the cosmological billairds on 
the UPHP ant with its restricted phase space is obtained by fixing a given energy at which the Hamiltonian flow is 
calculated. 

The operator te defined in Eq. EH can therefore be interpreted as the operator that, for each classical energy level 
E, acts on a semiclassical wavefunction (semiclassical in the sense that it is evaluated on a classical BKL configuration 
corresponding to a geiven sequence of epochs and eras). 

At each energy level E, the operator(s) te leave invariant (except for a ± sign) the eigenfunctions of the quantum 
eigenvalue problem; these Selberg eigenfunctions (??) are defined as generalizing the Riemann ^ function to closed 
(i.e. periodic) orbits (instead of prime numbers). 

The BKL map operators define the complete set of periodic orbits of cosmological billiards (from a classical point of 
view). 

The corresponding system of operators te can therefore be interpreted as the semiclassical operator which extracts 
an eigenvalue for each closed geodesics, which correspond to periodic orbit, according to it content of BKL epochs, 
BKL eras and the chosen symmetry-guotienting mechanism, and whose complete spectra are now classified. 


VIII. CONCLUDING REMARKS 

The aim of this investiagtion has been to define periodic phenomena for Cosmologial Billiards in 4 = 3 -I-1 spacetime 
dimensions. In particular, quantum BKL maps and the semiclassical BKL maps have eigenvalues which define the 
sign acquired by the wavefunctions for cosmological billairds, whose parity is defined by the corresponding periodicity 
phenomena. 

The paper is organized as follows. 

In the Introduction Cosmological Billiards are introduced. 

In Section IH the Billiard Modular Group is defined. 

The Big Billiard Group is defined in Section uni and the Small Billiard in Section CYl for which the differences with 
the Billiard Modular Group are outlined in Section El as those characterizing the symmetries of the solutions to the 
Einstein field equations. 

For this billiard groups, the language code is stated, and periodic orbits are spelled out and analyzed. 

The sign acquired by the Maass wavefunctions under the quantum BKL maps ?? characterizes the Selberg trace 
formula for Cosmological Billiards, which also defines the energy-level densities according to the quantum BKL maps 
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for periodic orbits, as analyzed in Scetion IVII As analyzed in Section IVII and initially commented in Section IVIIi 
the quantum BKL map operators are the same self-adjoint operator of the quantum BKL operators defined on the 
solution of the minisuperspace ruduction of the WDW equation for cosmological billiards. This result strictly applies 
to to the analysis of [i^ and and further specifies the analysis of and [i^ for Cosmological Billiards, 



FIG. 1. The domains of the big billiard V2{PGL{2,C)), the small billiard, the SL{2,Z) group, the SL{2,C) group and the 
billiard modular group. In particular, the domain of the big billiard is delimited by the geodesics u = 0, u = —1, + u — O', 

the domain of the small billiard, which coincides with that of the SL{2, Z) group, is delimited by the geodesics u = 0, w = —1/2 
and = 1; the domain of the SL{2,C) group is delimited by the geodesics u = 1/2, u = —1/2 and + v'^ = 1; the 

domain of the billiard modular group is delimited by the geodesics m = 0, u = l/2 and = 1. They are all plotted 

by solid black lines. The simmetry lines of the big billiard are represented by the dashed black lines. An oriented geodesic is 
drawn as gray dashed circle, and the oriented endpoints u'*' and u~ are indicated on the u axes. 
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